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The azimuthal collimation of di-hadrons with large rapidity separations in high multiplicity p+p 
collisions at the LHC is described in the Color Glass Condensate (CGC) effective theory [lj by N% 
suppressed multi-ladder QCD diagrams that are enhanced a^ 8 due to gluon saturation in hadron 
wavefunctions. We show that quantitative computations in the CGC framework are in good agree- 
ment with data from the CMS experiment on per trigger di-hadron yields and predict further 
systematics of these yields with varying trigger pr and charged hadron multiplicity. Radial flow 
generated by re-scattering is strongly limited by the structure of the p+p di-hadron correlations. In 
contrast, radial flow explains the systematics of identical measurements in heavy ion collisions. 



The discovery of di-hadron correlations in high multi- 
plicity proton-proton collisions [2J , long range in the an- 
gular (pseudo-rapidity) separation of the pairs relative to 
the beam axis and collimated in their relative azimuthal 
angle about this axis, provides significant insight into rare 
parton configurations in the proton and their dynamics 
in hadronic collisions. 

High multiplicity proton-proton collisions select "hot 
spot" configuations of wee gluon states in each proton. 
Quantum Chromodynamics (QCD) predicts that such 
hot spots have a maximum occupancy of order ag 1 
(as being the QCD fine structure constant), and have a 
typical size ~ 1/Qs, where Qs is a dynamical saturation 
scale. This scale grows with the energy and centrality of 
the collision; when Qs ^> Aqcd, the fundamental QCD 
scale, highly occupied hadron wavefunctions can be de- 
scribed using weak coupling methods. 

A weak coupling effective field theory (EFT) that de- 
scribes high density wee parton configurations in the pro- 
ton is the Color Glass Condensate (CGC) [5[. When 
CGC's shatter in a high multiplicity collision, multi- 
particle production is a consequence of approximately 
boost invariant radiation from "Glasma flux tubes" of 
transverse size l/Qs 0- Multiplicity distributions 
derived from factorization theorems [8j in this framework 
are in good agreement dEll with recent LHC data [ll|. 
Long range rapidity correlations of gluons computed in 
the CGC EFT [l2jwere previously shown to be in quali- 
tative agreement [l( with the CMS di-hadron correlation 
data. 

A source of long range rapidity correlations in hadron- 
hadron collisions are back-to-back gluons emitted from 
a single i-channcl gluon ladder; another source, called 
"Glasma graphs" are gluons emitted from two separate 
ladders. Representative graphs of each are shown in 
fig. (HJ). In the "dilute" high pr perturbative limit of 
QCD, the back-to-back contribution is dominant. How- 
ever, at high parton densities, when Q| >• Aq CD , and 
Pt ~ Qs> the effective coupling of gluons in ladders to 
strong color sources at higher rapidities changes from 



g — > 1/g. This corresponds to an enhancement of Glasma 
graphs by cv^ 8 compared to the a^ 4 enhancement of the 
back-to-back graphs. Equally important are the very 
different azimuthal structures of the two long range ra- 
pidity correlations. Back-to-back graphs, as the name 
suggests, are kinematically constrained to be peaked 
"away side" at relative azimuthal angle A(f> ~ it and have 
a negligible "near side" collimation at A(f> ~ as seen in 
high energy asym pto tics that produce a long range ra- 
pidity correlation 13|, [14 1 . 

In contrast, Glasma graphs give identical near and 
away side contributions because gluon emission is from 
independent ladders. As also noticed elsewhere IH ltjj |. 
these correlated contributions producing an azimuthal 
collimation are of order 1/N 2 ; their contribution would 
be negligible if one did not have the a^ 8 enhancement in 
the high multiplicity region. Within the CGC framework 
itself, there are leading N c correlated multi-ladder con- 
tributions [H which are long range in rapidity. However, 
these do not produce an azimuthal collimation 17, 18| . 
Likewise, there can be pomeron loop effects outside the 
framework of ref. Q ; again, these either do not give a col- 
limation or the collimation vanishes rapidly with increas- 
ing rapidity [l^jISOl- However, a Zero-Yield-at Minimum 



(ZYAM) procedure [21( used by the CMS collaboration 
to compute the per trigger near side di-hadron yield 22] 
only measures contributions that are collimated in Acp 
above the A0-independent background. The ZYAM pro- 
cedure allow one to focus on those di-hadron correla- 
tions that produced a collimated near side yield. This 
is also fortuitous because the relative normalization be- 
tween Glasma graphs and back-to-back graphs necessary 
to reproduce the di-hadron yield in the entire A$ range 
is not under theoretical control. 

The correlated two gluon glasma distribution, ex- 
pressed in terms of the two particle momentum space 



rapidities y p<q and transverse momenta p± , q± , is [12 1 
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FIG. 1. Representative back-to-back (left) and Glasma graphs 
(right) in perturbative QCD. 



where [23j Co 
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with D Am) = $ A2(1) (y„ q T + k T ) + $a 2(1) Qt - k T ). 
For our computation, we will also need the single inclu- 
sive gluon distribution 
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with the coefficient Ci 



4^(JVg-l)- 



The important in- 



gredient in these expressions is the unintcgrated gluon 
distribution (UGD) per unit transverse area, defined as 
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Here Ta is the forward scattering amplitude of a quark- 
antiquark dipole of transverse size rj_ on the target A; it, 
or equivalently, the UGD, is a universal quantity that 
can be determined by solving the Balitsky-Kovchegov 
(BK) equation [3, [IH as a function of the rapidity 
y = log(xo/a;). Solutions of the running coupling BK 
(rcBK) equation 26| used to compute structure func- 
tions are in good agreement with the HERA inclusive 
data (23|. 

The eqs. and (|3|) are obtained from the CGC for- 
malism in ref. §| for collisions of high high parton den- 
sity sources, as may be realized in nucleus-nucleus and 
high multiplicity proton-proton collisions. We emphasize 
that, albeit not shown explicitly in fig. (TT]), the deriva- 
tion of eq. (fT|) in ref. [l2| encodes the effect of radiation 
between the sources and the triggered gluons as well as 
the radiation between the gluons. In obtaining these re- 
sults, the distribution of color sources is assumed to be 
a non-local Gaussian distribution, whose variance is sim- 
ply related to ^a{v, k±). This assumed distribution has 
been shown recently to provide a good approximate solu- 
tion to the Balitsky-JIMWLK hierarchy for n-point light- 
like Wilson line correlators [2^, [2{|. The unintegrated 



gluon distribution in eq. (UJ have a "bell-shaped" struc- 
ture peaked at Qs, with the peak moving to larger k± 
with increasing rapidity. Thus eqs. (JI]) and ([3]) are in- 
frared finite unlike the expressions in ref. [7|. However, 
like the latter, they do not include multiple scattering 
contributions that are present in the formalism of ref. [8| 
and contribute for k± < Qs- Their effect on eq. (TTJ) is 
given by a non-perturbative constant [30j £ estimated nu- 
merically in ref. [3l| to be in the range 1/3-3/2. Fits 
to p+p multiplicity distributions for a range of energies 
at the LHC and A+A multiplicity distributions at RHIC 
obtained C = 1/6 E3- 

Given uncertainties in the 
numerical computation we will use the latter empirical 
value instead. 

Qualitatively, the origin of a large collimated contribu- 
tion from Glasma graphs occurs because the integrand in 
eq. ([T]) is large when the peaks of the "bell-shaped" un- 
integrated distributions strongly overlap, |kr| ~ |pr — 
kr| ~ |qT±kr| ~ Qs, giving a collimation at A$ = 0, tt. 
In practice, the result is smeared by fragmentation effects 
as well as the details of the integration. We will there- 
fore in the rest of this letter focus on the quantitative 
contribution of the Glasma graphs and compare the sys- 
tematics to the CMS proton-proton collision data. We 
will also explore the relative role of radial flow in generat- 
ing near side yields in proton-proton and nucleus-nucleus 
collisions. 

The initial condition in rcBK is the McLerran- 
Venugopalan-like (MV)- initial condition [32[ for Ta 
at the the initial rapidity (corresponding to Bjorken 
x = xq — 0.01). The minimum bias saturation scale Qq 
in the MV initial condition at the initial rapidity and the 
transverse area S± are adjusted to reproduce the single 
inclusive p+p distribution in eq. (]3])-for a more detailed 
discussion, see ref. Q . Diffractive scattering results from 
HERA indicate a strong dependence of the saturation 
scale on impact parameter or the centrality of the colli- 
sion [33|]. The centrality dependence of eq. ([JJ is there- 
fore studied here by keeping S± fixed and varying Qq at 
the initial rapidity scale 1341 . We follow the results of 
the HERA studies in ref. [33j and choose Qg values in 
the fundamental representation of 0.15, 0.3, 0.45 and 0.6 
GeV 2 to represent different centralities in computing the 
di-hadron yield [35| . 

Further calculational details are as follows. The strong 
coupling constant, as, is evaluated at the saturation 
scale. Because the di-hadrons of interest are widely sep- 
arated in rapidity, we assume that gluons fragment in- 
dependently. Fragmentation functions at forward rapidi- 
ties are not particularly well known (36j ; these will be 
better constrained as more data from the LHC becomes 
available. For our purposes, we consider a soft fragmen- 
tation function D\{x) = 3(1 — x) 2 /x and a hard one 
D2(x) = 2(1 — x)/x; the former is closer to the NLO fit 
function for gluon fragmentation to pions (37j , while the 
latter allows for hadrons to carry on average a larger frac- 
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FIG. 2. Associated yield for four different initial saturation 
scales representing different centralities. The blue circles (red 
squares) represent the softer (harder) D\(z) (D2(z)) fragmen- 
tation functions. Dashed lines connect computed points to 
guide the eye. Data points are from the CMS collaboration 



tion of the gluon momentum. Finally, we introduce an 
overall constant K factor which is the only parameter fit 
to the data presented here. It accounts for higher order 
corrections (both in our framework as well as in the frag- 
mentation functions) in addition to corrections necessary 
to fully implement the experimental acceptance [38| . 
For our analysis of the CMS data, we define 



trig 



dn i dV/^ dN fPT 



+2.4 <.p™ 

dr) 

-2.4 Jp™ in 







z 2 drjd 2 pT V z 



and 

d 2 N f +2A 

dK~4> =K 7- 2 .4 dr i? dr i* A (vp,v q ) 

r p ™" d P 2 T [ q ™ x dq 2 



(6) 



Ar\ and in the associated and trigger particle transverse 
momenta. As noted previously, the charged particle mul- 
tiplicity is varied by changing Qq in the initial conditions 
for rcBK evolution. We see that the agreement is quite 
good, especially at the highest multiplicities where we 
expect our formalism to perform best. At lower multi- 
plicities, the effect of high order corrections as well as 
impact parameter fluctuations become more important. 

In fig. ((3]), we plot the associated yield versus the pr 
trigger window for three associated particle windows as 
labeled. The top figure corresponding to 1.0 < p^ ssoc < 
2.0 GeV also shows the recent CMS measurements of 
the same quantity. Even though the overall normaliza- 
tion of our calculation needed to be augmented by a con- 
stant K-factor, K — 2.3, the momentum dependence of 
our results is parameter free. The other two plots are 
absolute predictions; though, as shown, they are quite 
sensitive to fragmentation. The sensitivity of the asso- 
ciated yield to different momentum cuts in our calcula- 
tion stems from an intrinsic scale (the saturation momen- 
tum) where the initial state wave-function is peaked. As 
argued in [l| maximal angular correlations occur when 
|Pt| ~ |qt| ~ Qs an d when pt and are parallel. This 
signal persists after including fragmentation and shows 
good agreement with the measured high multiplicity pp 
data. Our model computation provides strong support 
to the qualitative idea that a significant near side angu- 
lar correlation at semi-hard trigger and associate hadron 
momenta of 2— 4 GeV is evidence of saturation dynamics. 

Because the number of particles produced in the high- 
est multiplicity pp collisions are comparable to those in 
Cu Cu collisions one may speculate that (above and be- 
yond the collimation provided by our intrinsic QCD ef- 
fect) collective flow contributes significantly to the an- 
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employ a radial boost model where the angular distri- 
bution in A(j) in the laboratory frame is related to the 
corresponding distribution in A(f> in the local rest frame, 



TT minfmax 

Here p T 



-i min(max) 

and q T 



refer to bounds on 
the range of the trigger and associated hadron mo- 
menta respectively. Likewise, A?7 m i n (Arj max ) — 
2.0(4.0) denote the pseudo-rapidity gap [39( of hadrons 
within the experimental acceptance A{rj p ,rj q ) = 

(\Vp ~ Vq\ ~ ^Vrnin) (A?? max - \r} p - l] q \). 

The strength of the correlation in A(f> is quantified by 
the associated yield computed using the ZYAM proce- 
dure denned to be 
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where A0 m i n . is the angle at which the two particle cor- 
relation strength is minimal. In fig. ([2]), we plot the asso- 
ciated yield as a function of charged particle multiplicity, 
per minimum bias multiplicity, for the stated windows in 



where J is the Jacobian [42j relating distributions in the 
two frames. As transverse flow further collimates the 
signal, the overall strength of the associated yield will in- 
crease. However, the momentum dependence changes as 
well. The effect of the boost is demonstrated in fig. (j4}. 
Starting with our correlation in the local rest frame (red 
squares), we show the result after transverse boosts of 
(bottom to top) P = 0.1,0.2,0.25,0.3. One notices a 
qualitative change in the shape of the associated yield 
versus p^ lg - For smaller transverse boosts the depen- 
dence on Pp Ig is given by the intrinsic angular correlation 
generated by the Glasma graph of fig. . For large boost 
velocities the associated yield is driven by the collimation 
of the A<f> independent pedestal computed from the same 
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FIG. 3. Associated yield for central p+p (Q 2 = 0.6 GeV 2 ) col- 
lisions using soft (hard) Di (z)-blue circles (D2 (z)-red squares) 
fragmentation functions. Dashed lines connect computed 
points to guide the eye. The black squares are the vailable 
CMS data [HJ for the A > 110 multiplicity bin of pp col- 
lisions at ,/s = 7 TeV. The middle and bottom figures are 
predictions for the labeled associated pr windows. 



graph. Without any transverse flow this pedestal (or un- 
derlying event) is removed by the ZYAM procedure and 
therefore does not contribute to the associated yield. But 
after collimation, its signal exceeds that of the intrinsic 
angular correlation. The change in shape therefore sug- 
gests an upper bound of (3 = 0.25 in this simple model of 
flow in pp collisions. 

This is in complete contrast to heavy ion collisions 
where we expect flow to dominate the angular correla- 
tion [(| 0, 44 [. We demonstrate this with a compari- 
son of the p^ ls dependence of the (collimated by flow) 
pedestal in the Glasma with data from Pb Pb collisions 
at ^/s = 2.76 TeV in fig. ([5|). The agreement is quite 
good considering the very simple model of radial flow 
considered here. Flow effects here completely dwarf the 
intrinsic QCD correlations that were the dominant effect 
generating the near side azimuthal collimation in pp colli- 
sions. We should stress however that the pedestal (while 
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FIG. 4. Effect of transverse flow on the intrinsic pp correlation 
using the hard P>i{z) fragmentation function. Boosts from 
bottom to top: f3 = 0, 0.1, 0.2, 0.25, 0.3. 
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FIG. 5. Computations of the associated per trigger yield in 
Pb Pb collisions at yfs = 2.76 TeV using Q 2 , = 0.9 GeV 2 in the 
fundamental representation compared to the CMS data [2l| . 
The curves shown, with the P>i(z) fragmentation function, 
are for transverse boosts of ft = 0, 0.65, 0.85. At large flow 
velocities, the intrinsic angular correlation is entirely washed 
out. 



independent of A<fi in the local rest frame) is also an in- 
trinsic two-particle correlation generated by the Glasma 
graph and the p^ lg dependence seen in fig. ([5]) is repre- 
sentative of this underlying dynamics. 
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